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摘 要 : 地 表 粗 糙 度 作为 风蚀 防护 中 常用 的 传统 指标 ,在 防护 林 、 沙 障 等 工程 措施 的 防护 效益 评价 


中 发 挥 了 重要 作用 。 但 在 实际 应 用 中 发 现 ,地 表 粗 糙 度 存在 不 完全 准确 的 现象 。 针 对 这 一 问题 ， 


引入 新 的 参考 指标 附 面 层 位 移 厚度 ,借助 风 洞 模拟 实验 的 方法 ,通过 测定 不 同 讽 长 聚 乳 酸 沙袋 沙 


障 中 的 风速 频数 、 拟 合 风速 流 场 . 计 算 防 风 效 能 等 
的 准确 性 。 结 果 表 明 : 随 着 聚 乳酸 沙袋 沙 障 障 格 边 长 的 增 大 ,平均 风速 呈现 出 上 升 的 趋势 ,其 中 障 
格 讽 长 为 1.$m 和 2.0 冲 沙 障 内 的 平均 风速 分 别 是 1.0 征 沙 障 的 1.13 和 1.24 倍 ;防风 效能 则 呈现 出 
降 的 趋势 ,其 中 障 格 边 长 为 1.0 m 沙 障 内 各 处 均 大 于 0.6,1.5 m 沙 障 内 各 处 均 大 于 0.5,2.0 m^ 
各 处 均 大 于 0.4; 风 速 流 场 高 速 区 的 面积 呈现 出 增 大 的 趋势 。 这 表明 随 着 沙 障 障 格 讽 长 的 增 大 ， 
沙 障 的 防护 作用 逐渐 下 降 。 与 此 同时 , 随 着 沙 障 障 格 边 长 的 增 大 , 附 面 层 位 移 厚度 呈现 出 减 小 的 


a 二 


来 验证 地 表 粗 糙 度 和 附 面 层 位 移 厚 度 这 2 个 指标 


趋势 ,而 地 表 粗 糙 度 则 呈现 出 先 减 小 后 增 大 的 趋势 。 综 上 , 附 面 层 位 移 厚 度 在 评价 沙 障 防护 作用 


时 的 准确 性 更 高 ,可 作为 新 的 补充 参考 指标 。 


关 d$ 词 : 附 面 层 位 移 厚度 ; 风蚀 评价 ; 风 洞 模拟 ; 聚 乳酸 沙 障 
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风力 侵蚀 是 指 在 气流 冲击 作用 下 土 粒 、 沙 粒 脐 
离 地 表 ,被 搬运 和 堆积 的 过 程 "。 风 蚀 造 成 了 严重 
的 自然 灾害 和 环境 问题 ,尤其 是 在 干旱 半 干 旱地 
区 ,往往 成 为 制约 农业 生产 和 导致 环境 恶化 的 主导 
因素 ,影响 了 这 些 地 区 的 资源 开发 和 社会 经 济 的 持 
续 发 展 ” 。 为 有 效 降低 风沙 灾害 ,研究 人 员 开展 了 
多 种 风蚀 防护 措施 的 研究 ,如 防护 林 、 生 物 结 皮 、 固 
沙 剂 及 沙 障 等 4“”。 其 中 , 沙 障 因 其 布设 简便 、 防 护 
效果 好 及 防护 时 间 长 等 优点 ,是 我 国防 沙 治 沙 工作 
中 的 重要 工程 措施 ”“”。 在 风蚀 防护 措施 的 防护 效 
益 评 价 中 ,地 表 粗 糙 度 因 其 对 下 垫 面 改变 的 敏感 度 
最 高 ,得 到 了 广泛 的 应 用 """。 
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由 于 地 表 起 伏 以 及 粗糙 不 光滑 的 性 质 ,气流 流 
经 地 表 时 受 摩 擦 力 的 影响 ,导致 近 地 表 处 风速 降 
低 ,在 垂直 高 度 上 呈现 出 不 同 分 布 规律 ,这 种 风速 
随 高 度 分 布 的 曲线 称 为 风速 廓 线 “。 在 风速 廓 线 
上 风速 为 零 的 位 置 并 不 出 现在 高 度 为 零 的 地 表 ,而 
是 高 出 一 定 距 离 ,这 一 高 度 被 称 为 地 表 粗 糙 度 ,又 
称 空气 动力 学 地 表 粗 烽 度 。 近 代 流 体力 学 奠基 
人 普 衣 特 指出 粗糙 度 概念 的 本 质 表 明了 气体 具有 
粘性 ,并 给 出 了 粗糙 地 面 的 关系 式 ,后 经 尼古拉斯 
的 实验 测定 ,得 出 普 适 速度 分 布 公式 “: 


hV, 
u= v(sas wess) (1) 
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IIP iu WIRE (mss); ,为 摩 阻 流速 (ms D 5h DJ 
BJE Cm) vis SIR E CO s) 

随后 风沙 物理 学 黄 基 人 和 拜 格 诺 将 上 述 理论 引 
入 到 风沙 运动 中 ,并 提出 了 固定 地 表 的 风速 廓 线 
go 


u, =5.75V, IE 元 (2) 
0 


式 中 : u, 为 高 度 h 处 的 风速 (ms 1); Z, 为 粗糙 度 (m ) 
进一步 的 ,可 以 推导 出 粗糙 度 的 计算 公式 : 


272 
lgZ,— 7 182 
pee t (3) 


AP: Z,.Z, 为 2 个 不 同 的 高 度 (m); v, V v, 分 别 为 
高 度 Zi 和 2, 处 的 风速 (ms')。 

地 表 粗 糙 度 反映 了 粗糙 地 表 对 气流 的 阻 涉 作 
用 ,表明 地 表 不 同 下 垫 面 与 近 地 表 气流 的 相互 作用 
和 物质 能 量 交 换 , 在 风蚀 防护 评价 中 得 到 了 广泛 的 
应 用 "“* i。 但 在 近 几 十 年 的 国内 外 研究 中 ,科研 人 
员 大 多 只 是 采纳 这 一 参数 ,并 没有 进行 更 深入 的 研 
究 , 甚 至 有 学 者 提出 现 有 的 计算 方式 只 是 对 粗糙 度 
的 估算 ,计算 结果 并 不 准确 后。 通过 分 析 粗 糙 度 的 
计算 公式 ,本 研究 发 现 ,粗糙 度 由 2 个 不 同 高 度 的 风 
速 计算 获 得 ,而 不 同 的 高 度 组 合 会 得 到 不 同 的 数 
值 ,因此 粗糙 度数 值 具 有 不 唯一 性 ,进而 只 适用 于 
当 次 实验 ,无 法 用 于 多 实验 间 的 对 比分 析 。 地 表 粗 
烟 度 在 特定 下 垫 面条 件 下 并 不 是 一 个 常数 ,而 是 随 
着 风速 的 变化 产生 明显 的 变化 ,野外 气流 变化 较为 
复杂 , 仅 用 2 个 高 度 的 数据 进行 计算 所 得 到 的 误差 
也 较 大 。 

通过 对 风速 廓 线 的 研究 ,本 研究 提出 利用 附 面 
层 位 移 厚 度 来 进行 评价 防护 效益 。 附 面 层 是 指 由 
于 气流 具有 粘性 而 不 可 忽略 的 地 表 流 动 层 ™, 附 面 
层 位 移 厚 度 是 指 相 比 于 光滑 通道 ,相同 体积 的 粘性 
流体 在 相同 时 间 内 流 过 的 通道 所 需 加 宽 的 厚度 ”2 ， 
代表 了 粗糙 地 表 对 气流 能 量 的 消耗 。 将 风速 廓 线 
进行 积分 ,计算 函数 曲线 与 坐标 轴 间 的 面积 ,与 其 
面积 相等 的 矩形 的 长 度 即 为 附 面 层 位 移 厚度 (图 
1), 进 一 步 可 推导 出 附 面 层 位 移 厚 度 的 计算 公式 

后 we 

式 中 : Ad 为 附 面 层 位 移 厚 度 ; f(x) 29 DAUR ERA PRI 
数 ; a 为 理想 状态 下 的 风速 值 (ms ')。 


(4) 


Ad 高 度 /m 


1 附 面 层 位 移 厚 度 示意 图 
Fig. 1 Diagram of displacement thickness of boundary layer 


聚 乳酸 沙袋 沙 障 是 近年 来 应 用 较 多 的 新 材料 
沙 障 , 具 有 绿色 可 降解 .防护 时 间 长 成 本 低 ,布设 
简单 等 优点 , 且 实现 了 以 沙 治 沙 的 理念 ,在 风沙 治 
理 中 发 挥 着 重要 的 作用 一 ”。 因 此 ,本 研究 以 聚 乳 
酸 沙袋 沙 障 ( 图 2) 为 研究 对 象 来 分 析 地 表 粗 糙 度 和 
附 面 层 位 移 厚 度 在 评价 防护 措施 防护 效果 时 的 准 
确 性 ,通过 风 洞 模拟 实验 ,分别 计算 地 表 粗 糙 度 和 
附 面 层 位 移 厚 度 ,并 以 风速 频数 .风速 流 场 及 防风 
效能 等 指标 为 依据 来 衡量 这 2 个 参数 的 准确 性 和 适 
用 性 。 


1 材料 与 方法 


1.1 RUA 
本 研究 中 应 用 的 风 洞 模拟 实验 系统 属于 直流 
式 风 洞 ( 图 3a) 。 风 洞 全 长 约 为 24.3 m ,其 中 实验 段 
长 度 为 12 m, 实 验 段 的 横 截 面积 为 0.6 mx0.6 m; 动 
力 段 长 度 为 1.48 m ,为 实验 提供 气流 ;稳定 段 长 度 为 
1.5 m, 内 置 蜂 窜 网 ,将 大 团 空气 打 散 ,转变 为 细 流 ; 
收缩 段 长 度 为 1.5 m, 通 过 横 截面 积 的 急剧 收缩 加 速 
气流 速度 ,三 者 可 提供 风速 在 3~40 ms 范围 内 
连续 可 调 的 稳定 气流 。 在 实验 中 ,我 们 利用 粗糙 元 
(4 emX5 cmx6 cm, 图 3b) 来 拟 合 风速 廊 线 ,布设 方 
式 如 图 3c 所 示 , 并 通过 以 下 公式 计算 雷诺 数 (Re)™. 
Re= 4*4 (5) 
式 中 :wu 为 流速 (ms ');d 为 特征 长 度 (m) ;vw 为 空气 
粘度 (ms ), 取 值 15x10? ms 。 人 参考 野外 气象 
资料 ,本 研究 设置 实验 风速 为 8 ms CE 1 mAb), 
此 时 雷诺 数 为 3.2x10’。 
为 准确 测量 不 同位 置 的 风速 ,我 们 采用 了 三 维 
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以 聚 乳 酸 (PLA ) 沙 障 为 例 1975 


图 2 聚 乳酸 沙袋 沙 障 
Fig. 2 Polylactic acid sandbag barrier 


移 测 系统 ,可 在 实验 段 自由 移动 ,精度 为 1 mm。 测 
量 设备 我 们 采用 了 热 膜 风速 仪 (IFA300,TSI, Shorev- 
iew, MN, USA ) 和 热线 风速 仪 (VT-200,KIMO , Mont- 
pon Ménestérol, France) ; 
1.2 模型 制作 

聚 乳 酸 沙 袋 沙 障 常用 障 体高 度 为 10 cm, 障 格 
边 长 为 1.0m、1.5m 和 2.0m。 本 研究 采用 聚 乳酸 材 
料 灌 装 沙 土 的 方式 按照 1:5 的 比例 制作 模型 ,模型 
高 度 为 2 em, 相应 障 格 边 长 分 别 为 20 cm ,30 cm 和 
40 cm , 铺 成 10 行 3 列 网 格 ,此 外 为 接近 野外 地 表 实 
际 状况 ,在 实验 段 底板 上 粘 附 一 层 沙 粒 。 

风 洞 实验 测 点 布设 如 图 4 所 示 。 沿 风向 在 洞 体 
中 间 网 格 的 中 线 布设 测 点 ,初始 位 置 位 于 障 前 20 em, 
水 平方 向 测 点 间隔 为 障 格 的 十 等 分 中 ,垂直 高 度 测 
点 间隔 为 1 em, 垂直 高 度 测量 范围 为 1~20 cmo TI 
风向 测定 不 同 障 格 内 风速 ,直至 网 格 内 风速 值 趋 于 
稳定 ,并 在 下 一 个 网 格 内 布设 水 平 截面 测 点 。 水 平 
截面 测 点 每 行 . 列 分 别 布设 9 个 测 点 ,布设 高 度 分 别 
为 1 cm( 低 于 障 高 )、2 cm( 与 障 高 平 齐 ) 和 4 cem 
于 障 高 )。 
1.3 风速 频数 统计 

统计 风速 稳定 后 障 格 内 不 同 高 度 的 风速 频数 ， 
分 析 最 小 值 . 最 大 值 平均 值 .标准 差 . 峰 度 isi HE e 
变异 系数 等 指标 。 峰 度 是 表征 与 正 态 分 布 曲线 相 
比 ,概率 密度 分 布 曲线 分 布 的 陡峭 程度 , 正 值 表示 
平均 值 处 分 布 更 多 的 数值 , 负 值 表示 数值 分 布 较为 
集中 ,计算 公式 为 


1S\@,-2) 
KS (6) 
Ca 
式 中 :有 为 峰 度 ; x, 为 第 ;个 统计 数值 ; x 为 样本 平均 
值 ;为 样本 数量 。 

偏 度 是 指 相对 于 正 态 分 布 曲线 ,概率 密度 分 布 
曲线 分 布 的 对 称 情况 。 偏 度 为 正 值 表明 数据 更 多 
的 分 布 在 左 侧 ,样本 中 位 数 小 于 均值 ; 偏 度 为 负 值 
表明 数据 更 多 的 分 布 在 右 侧 ,样本 中 位 数 大 于 均 
值 。 计 算 公式 为 


Be (7) 


式 中 :5 为 偏 度 ;x 为 第 i 个 统计 数值 ; x 为 样本 平均 
值 ;n 为 样本 数量 。 

样本 峰 度 和 侦 度 可 以 分 析 样 本 频数 分 布 特征 ， 
变异 系数 则 可 以 分 析 样 本 的 离散 程度 。 一 般 认为 ， 
数值 大 于 0.15 时 样本 数据 变异 程度 较 大 ,数据 可 徘 
性 较 差 。 计 算 公 式 为 ”， 


Er (8) 
式 中 : C, 为 变异 系数 ;sd 为 标准 差 ; mean 为 样本 


均值 。 
1.4 空间 相关 性 

空间 自 相关 统计 量 是 用 于 度量 空间 数据 的 一 
个 基本 性 质 , 即 某 位 置 上 的 数据 与 其 他 位 置 上 的 数 
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(a) 洞 体 示意 图 
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(b) 粗糙 元 示意 图 
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(c) 粗糙 元 摆 放 位 置 
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图 3 风 洞 示意 图 
Fig. 3 Wind tunnel schematic diagram 


据 间 的 相互 依赖 程度 ,模型 包括 块 金 值 . 基 台 值 、 变 
程 等 参数 。 其 中 , 变 程 指 相 邻 2 个 位 置 具 有 空间 相关 
性 的 最 大 距离 , 块 金 值 与 基 台 值 的 比值 即 为 空间 相 
KJE ,一般 认 为 区 域 中 的 测 点 间 空 间 相关 度 在 0~0.25 
范围 内 具有 极 强 的 相关 性 ,在 0.25~0.75 范 围 内 具有 
中 等 程度 的 相关 性 ,大 于 0.75 时 空间 相关 性 较 弱 *。 
1.5 防风 效能 

防风 效能 是 指 气流 经 过 防护 措施 后 风速 下 降 
的 程度 。 计 算 公式 为 2 : 


gru i 


(9) 
Ug 
式 中 :a 为 防风 效能 ; w AKU Gms!) ; v, 为 气流 
经 过 防护 措施 后 的 风速 (ms )。 
2 结果 与 分 析 


21 地 表 粗 糙 度 与 附 面 层 位 移 厚 度 
通过 测量 风速 达到 稳定 后 的 障 格 内 不 同 高 
度 处 的 风速 值 , 拟 合 出 不 同 边 长 沙 障 的 风速 廓 线 
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以 聚 乳 酸 (PLA ) 沙 障 为 例 1977 
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图 4 实验 测 点 分 布 


Fig. 4 Distribution of experimental measurement points 


(图 5), 各 风速 廊 线 均 能 较 好 地 拟 合 为 以 自然 指数 
为 底 的 指数 函数 ,上 且 决定 系数 ( 尼 ) 均 大 于 0.97。 

根据 风速 廊 线 可 以 计算 出 不 同 边 长 的 障 格 内 
风速 稳定 后 的 粗糙 度 和 附 面 层 位 移 厚度 ( 表 1) ,其 
中 在 地 表 粗 糙 度 的 计算 中 选取 了 实际 高 度 为 1$ cm 
和 100 em 处 的 风速 值 。 从 表 1 可 见 , 地 表 粗 糙 度 随 
着 障 格 边 长 的 增加 呈现 出 先 减 小 后 增 大 的 趋势 , 边 
长 为 1.0mm 沙 障 网 格 的 地 表 粗 烙 度 数值 最 高 ,分 别 
为 15Sm 和 2.0m 的 2.98 和 1.80 倍 。 与 地 表 粗 糙 度 不 
同 的 是 , 附 面 层 位 移 厚 度 呈 现 出 随 着 障 格 边 长 的 增 
大 而 减 小 的 规律 。 
2.2 风速 频数 统计 

统计 分 析 不 同 边 长 沙 障 风速 稳定 后 的 障 格 内 
风速 频数 ( 表 2) 可 见 ,不 同 高 度 的 风速 频数 均值 随 
着 障 格 边 长 的 增 大 均 呈 现 出 上 升 的 趋势 。 边 长 为 
1.5 m 4112.0 m 的 沙 障 障 格 内 风速 频数 均值 在 低 于 障 
格 高 度 处 分 别 是 1.0 m 障 格 的 1.28 和 1.69 倍 ,在 障 
格 高 度 处 分 别 是 1.0 m 障 格 的 1.14 和 1.20 倍 ,在 高 
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于 障 格 高 度 处 分 别 是 1.0 m BER AY 1.06 RI 1.114%, 
整体 平均 分 别 是 1.0m 障 格 的 1.13 和 1.24 倍 。 从 增 
长 的 幅度 可 以 看 出 , 随 着 高 度 的 增加 ,风速 趋向 于 
稳定 。 但 是 数据 标准 差 较 大 ,表明 均值 受 极 值 的 影 
响 较 大 。 且 从 数据 分 布 的 峰 度 、 偏 度 可 以 看 出 , 风 
速 频数 离散 程度 较 大 ,表明 障 格 内 的 风速 变化 较为 
UU HER. HEP , 近 半 数 样本 的 变异 系数 达到 
0.15 以 上 ,这 表明 风速 均值 受 极 值 影响 较 大 ,需要 更 
进一步 的 处 理 和 分 析 。 
2.3 空间 自 相关 性 

为 了 进一步 研究 沙 障 防 护 效果 和 障 格 边 长 之 
间 的 关系 ,需要 拟 合 障 格 内 的 风速 流 场 。 在 拟 合 之 
前 ,为 分 析 障 格 内 测 点 之 间 是 否 具 有 空间 相关 性 ， 
利用 地 统计 学 的 方法 对 不 同 边 长 的 障 格 风 速 频数 
进行 空间 自 相 关 分 析 ( 表 3)。 结 果 表 明 ,不同 边 长 
的 沙 障 网 格 内 的 风速 频数 均 能 拟 合 为 高 斯 模型 H. 
空间 相关 度 均 小 于 0.25 ,表明 障 内 测 点 间 具 有 极 强 
的 空间 相关 性 , 风 洞 模拟 实验 数据 具有 较 高 的 可 靠 
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5 风速 稳定 后 沙 障 网 格 内 风速 廓 线 
Fig. 5 Wind speed profile within the sand barrier grid after stable wind speed 
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表 1 不 同 障 格 边 长 的 沙 障 地 表 粗 糙 度 和 
附 面 层 位 移 厚 度 
Tab.1 Surface roughness and boundary layer 


displacement thickness of sand barriers with 


different side lengths 
参数 障 格 边 长 
1.0 m 1.5m 2.0 m 
地 表 粗 糙 度 /m 0.0137 0.0046 0.0076 
附 面 层 位 移 厚度 20.00 19.38 18.09 
性 和 准确 性 。 


2.4 风速 流 场 

数值 模拟 不 同 障 格 内 风速 稳定 后 流 场 如 图 6 所 
示 。 不 同 边 长 障 格 内 风速 均 小 于 原 风 速 ,表明 沙 障 
起 到 了 良好 的 防护 作用 。 测 量 高 度 为 5 em 处 的 风 
速 流 场 可 以 看 到 , 边 长 为 1.0 m 的 障 格 在 障 后 左右 
形成 2 个 蓝 色 减速 区 ,中 间 形 成 黄色 加 速 区 ,核心 有 
部 分 橙色 区 域 ; 边 长 为 1.5m 障 格 后 形成 较 大 范围 
的 蓝 色 低速 区 ,随后 风速 逐渐 恢复 ,并 在 障 格 后 半 


uq E 
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TAE IX, 2 PE CEH D TEE DX Sa, ; 边 长 为 2.0m 障 
格 在 障 后 1/3 范围 内 形成 2 个 蓝 色 低 速 区 ,中 后 部 则 
是 大 范围 的 橙色 加 速 区 。 测 量 高 度 为 10 em 的 风速 
流 场 可 以 看 到 , 边 长 为 1.0m 惠 障 格 后 风速 逐渐 降低 ， 
两 侧 速 度 降低 更 快 ; 边 长 为 1.5 m 障 格 内 形成 中 部 
低速 ,向 两 侧 增 大 的 结构 ;2.0 m 障 格 内 流 场 结构 与 
测量 高 度 5 em 时 相似 ,但 颜色 变 浅 。 测 量 高 度 为 
20 cm 处 不 同 边 长 障 格 均 呈现 出 障 后 风速 下 降 的 趋 
势 ,日 障 格 边 长 越 大 ,风速 越 大 。 
2.5 防风 效能 

通过 分 析 风 速 稳定 后 不 同 边 长 沙 障 障 格 内 的 
防风 效能 (图 7) 可 见 , 不 同 边 长 障 格 内 的 防风 效能 
均 大 于 0.4, 且 随 着 边 长 的 增 大 ,防风 效能 呈现 出 下 
降 的 趋势 。 边 长 为 1.0 m 障 格 内 的 防风 效能 均 大 于 
0.6, 边 长 为 1.5 m 障 格 内 的 防风 效能 均 大 于 0.5, 边 
长 为 2.0 m 障 格 内 防风 效能 均 大 于 0.4。 边 长 为 1.0m 
和 1.5 m 的 障 格 内 防风 效能 下 降 较为 均匀 , 边 长 为 
2.0 m 的 障 格 内 防风 效能 在 0.4~0.6 范 围 内 急剧 下 


表 2 不 同 边 长 沙 障 障 格 内 风速 频数 分 布 


Tab. 2 Frequency distribution of wind speed within different side length barriers 


障 格 边 长 /m 高 度 /em 最 小 值 /ms 最 大 值 /ms 均值 + 标准 差 /m.s” IE BE DIES 变异 系数 
1.0 2 0.69 3.02 1.86x0.62 -1.07 —-0.12 0.33 
10 2.55 4.93 3.62x0.55 -0.42 0.17 0.15 
20 4.01 6.75 4.89+0.42 3.45 0.96 0.09 
1.5 5 0.46 3.99 2.39+1.05 -].25 —-0.31 0.44 
10 3.05 5.60 4.13+0.52 0.02 0.24 0.13 
20 4.36 5.98 5.20+0.42 0.12 0.37 0.07 
2.0 5 0.70 4.55 3.15+1.05 -0.21 —0.98 0.33 
10 2.97 4.93 4.33+0.41 0.56 -0.72 0.10 
20 4.53 6.88 5.42+0.35 —0.68 0.15 0.08 
表 3 不 同 边 长 的 沙 障 障 格 空间 自 相 关 分 析 
Tab.3 Autocorrelation analysis of barrier grid space with different side lengths 
障 格 边 长 /m 模型 块 金 值 基 台 值 变 程 空间 相关 度 决定 系数 (RR) 残 差 平方 和 (RSS)/107 
1.0 Gaussion 0.045 0.505 6.85 0.09 0.998 0.751 
Gaussion 0.064 0.371 7.43 0.17 0.998 0.603 
Gaussion 0.056 0.285 11.22 0.20 0.993 0.695 
1.5 Gaussion 0.048 1.071 8.80 0.04 0.999 1.181 
Gaussion 0.035 0.369 TAT 0.09 0.998 1.115 
Gaussion 0.001 0.130 3.37 0.01 0.782 7427 
2.0 Gaussion 0.047 0.436 4.55 0.01 0.958 7.330 
Gaussion 0.036 0.374 5.67 0.10 0.997 4.395 
Gaussion 0.057 0.673 7.23 0.08 0.991 7.691 
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(a) 边 长 1.0 m 障 格 , 测量 高 度 5 cm (b) 边 长 1.0 m 障 格 , 测量 高 度 10 cm (c) 边 长 1.0 m 障 格 , 测量 高 度 20 cm 
20 20 
18 18 
16 16 
14 14 
g g 12 g 12 
o o o 
id gj 10 ij 10 
a ii 8 dg 
6 6 
4 Si) 4 
2 BN A 2 
0 0 
2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 20 
位 置 /cm 位 置 /cm 位 置 /cm 
(d) 边 长 1.5 m 障 格 , 测量 高 度 5 cm (e) 边 长 1.5 m 障 格 , 测量 高 度 10 cm (f) 边 长 1.5 m 障 格 , 测量 高 度 20 cm 
30 
25 
20 
E E E 
机 Em "BE 
ii al ii 
10 
5 
0 
5 10 15 20 25 30 5 10 15 20 25 30 5 10 15 20 25 30 
位 置 /cm 位 置 /cm 位 置 /cm 
(g) 边 长 2.0 m 障 格 , 测量 高 度 5 cm (h) 边 长 2.0 m 障 格 , 测量 高 度 10 cm (i) 边 长 2.0 m 障 格 , 测量 高 度 20 cm 
0 
5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 
位 置 /cm 位 置 /cm 位 置 /cm 
en NSPT 
o x 9 «& $6 $ x o d 6 9 «x € d 6 o « o 
O Oc CO — TFN c c e e sr TFT FH ww ND CO ND 


图 6 不 同 边 长 障 格 风速 流 场 
Fig.6 Wind speed flow fields of different side length barrier grids 


降 ,表明 边 长 较 大 的 网 格 防风 蚀 能 力 较 差 ,网 格 内 
风速 较 高 。 
2.6 相关 性 分 析 

地 表 粗 糙 度 和 附 面 层 位 移 厚 度 与 各 指标 间 的 
相关 性 分 析 ( 图 8) 可 见 , 地 表 粗 糙 度 与 障 格 边 长 和 
风速 呈 负 相关 关系 ,与 防风 效能 呈正 相关 关系 ,但 
相关 性 不 显著 。 附 面 层 位 移 厚 度 与 障 格 边 长 和 风 
速 呈 负 相 关 关 系 , 与 防风 效能 呈正 相关 关系 , 且 相 


关 性 极 显 车 ,相关 系数 均 大 于 0.96。 
3 讨论 


3.1 地 表 粗 糙 度 的 准确 性 

由 表 2、 图 6 和 图 7 可 见 , 聚 乳酸 沙袋 沙 障 的 防 
护 效 果 随 着 障 格 边 长 的 增 大 而 减弱 。 从 表 1 可 以 看 
到 , 附 面 层 位 移 厚度 这 一 指标 准确 地 反映 出 沙 障 防 
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图 7 不 同 边 长 沙 障 防风 效能 
Fig. 7 Wind prevention efficiency of sand barriers with 


different side lengths 


护 效 果 的 规律 性 ,地表 粗糙 度 的 结果 并 不 完全 准 
确 , 且 图 8 相关 性 分 析 也 表明 地 表 粗 糙 度 与 各 参数 
间 的 相关 性 低 于 附 面 层 位 移 厚 度 。 为 了 验证 地 表 
粗糙 度 的 适用 性 ,本 研究 选取 多 个 高 度 的 风速 值 与 
实际 1m 高 度 处 的 风速 组 合计 算 地 表 粗 糙 度 ( 表 
4)。 从 表 4 可 见 ,实际 计算 高 度 在 15~20 cm 时 ,地 
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图 8 各 指标 间 相 关 性 分 析 


Fig. 8 Correlation analysis between various indicators 


表 粗 糙 度 并 没有 反应 出 真实 情况 ,而 其 他 高 度 计算 
出 的 结果 则 有 较 好 的 指示 性 。 

本 研究 进一步 拟 合 出 沙 障 障 格 沿 风向 的 纵 截 
面 流 场 图 (图 9)。 由 图 9 可 见 , 气 流 在 经 过 沙 障 时 出 
现 明 显 的 抬升 加 速 作 用 , 沙 障 上 方 区 域内 抬升 气流 
与 原 空气 看 加 作用 明显 ,风速 变化 较为 剧烈 ,而 在 


表 4 不 同 高 度 组 合计 算 所 得 的 地 表 粗 糙 度 


Tab. 4 Surface roughness calculated by combining different heights 


nea Pm 5cm 10 cm 15 cm 
1.0 0.0198 0.0243 
1.5 0.0092 0.0094 
2.0 0.0022 0.0059 


注 :不 同 小 写字 母 表示 在 P=0.05 水 平 上 差异 显著 。 
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图 9 沿 风 癌 纵 截面 流 场 


Fig.9 Flow fields of longitudinal section along the wind direction 
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抬升 区 域 之 上 ,气流 又 趋 于 稳定 。 表 4 中 地 表 粗 糙 
度 异 常 值 的 计算 高 度 正好 对 应 于 模型 中 气流 急剧 
变化 的 3~4 cm 的 高 度 , 因 此 可 以 判定 气流 的 剧烈 变 
化 是 导致 地 表 粗 糙 度 结果 不 准确 的 主要 原因 。 
3.2 附 面 层 位 移 厚 度 的 适用 性 

在 实际 中 ,地表 附 面 层 的 高 度 要 远 超 出 防护 措 
施 的 有 效 防 护 高 度 , 附 面 层 位 移 厚度 也 超过 本 研究 
的 测量 范围 ,因此 在 应 用 于 风蚀 防护 评价 时 需要 进 
行 一 定 的 改进 。 在 式 (4) 中 ,理想 状态 下 的 风速 值 
(ao) 无 法 测量 ,而 在 应 用 中 选取 超过 防护 措施 有 效 
防护 高 度 的 风速 值 即 可 。 由 图 9 可 见 ,不同 边 长 的 
沙 障 障 格 内 实际 高 度 超过 0.5 m 时 ,风速 趋 近 于 原 
风速 ,因此 可 以 选取 1.0 m 作 为 评价 沙 障 防护 效果 
时 附 面 层 位 移 厚 度 的 计算 高 度 。 此 外 ,在 进行 防护 
林 、 生 物 结 皮 等 其 他 防护 措施 的 效果 评价 时 也 可 以 
进行 相应 调整 。 


4 结论 


(1) 聚 乳 酸 沙袋 沙 障 具有 较 好 的 防风 作用 , 且 
防护 效果 随 着 沙 障 规 格 的 增 大 而 减弱 。 

(2) 地 表 粗 糙 度 是 由 2 个 不 同 高 度 的 风速 计算 
得 到 , 当 高 度 选取 在 风速 急剧 变化 的 区 间 范 围 内 
时 ,计算 得 到 的 地 表 粗 糙 度 准确 性 较 差 ,因此 在 应 
用 时 要 注意 气流 结构 特征 。 

(3) 附 面 层 位 移 厚度 是 对 有 效 防护 高 度 内 的 风 
速 廊 线 进行 计算 ,弥补 了 区 间 气 流 急剧 变化 带 来 的 
影响 ,在 评价 防护 措施 的 效果 时 准确 性 较 高 ,可 以 
作为 地 表 粗 糙 度 的 辅助 参考 指标 。 
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Applicability study of the boundary layer displacement thickness to the 
evaluation of sand barrier windbreak benefits: A case of 
polylactic acid sand barrier 
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Abstract: Surface roughness, a conventional index widely applied in wind erosion protection, plays a crucial role 
in evaluating the effectiveness of engineering measures such as forest and sand barriers. Despite its prevalent use, 
practical applications reveal its limited accuracy. Consequently, we introduce a novel index termed “boundary lay- 
er displacement thickness”. Through wind tunnel simulations, we assessed the wind speed frequency distribution 
of various sandbag barriers. We computed the mean wind speed, wind prevention efficiency, and also scrutinized 
the wind speed flow field to validate the precision of both roughness and boundary layer displacement thickness 
measurements. The findings indicate that an increase in the length of the polylactic acid sandbag barrier corre- 
sponds to an upward trend in average wind speed. Specifically, the average wind speed within the 1.5-m and 2.0-m 
barriers is 1.13 and 1.24 times that of the 1.0-m barrier, respectively. Concurrently, wind prevention efficiency ex- 
periences a decline, with values exceeding 0.6 for the 1.0-m barrier, 0.5 for the 1.5-m barrier, and 0.4 for the 2.0-m 
barrier. Analysis of the wind velocity field data reveal that the high-speed zone expands as the barrier grid size in- 
creases, signifying a weakened protective effect with large sand barriers. Notably, the boundary layer displace- 
ment thickness demonstrates a consistent decreasing trend with increasing barrier size, aligning with the afore- 
mentioned pattern. However, surface roughness exhibits an anomalous pattern, initially decreasing and then in- 
creasing. In summary, the boundary layer displacement thickness proved to be more accurate and emerged as a 
promising supplementary reference indicator for evaluating the efficacy of wind erosion protection measures. 

Key words: boundary layer displacement thickness; wind erosion evaluation; wind tunnel simulation; polylac- 


tic acid sand barrier 


